Preprint DFPD 02/TH/04 (Univ. of Padua) 1-26 (2008) 



Three-Body Dynamics in One Dimension: 
A Test Model for the 3-Nucleon System 
with Irreducible Pionic Diagrams 

T. Melde \ L. Canton 1 ' 2 , J.P. Svenne 3 



1 Dipartimento di Fisica, Universita di Padova, Via F. Marzolo 8, Padova, Italy, 1-35131 

2 Istituto Nazionale di Fisica Nucleare, Sezione di Padova, Via F. Marzolo 8, Padova, Italy, 
1-35131 

3 Winnipeg Institute for Theoretical Physics and Department of Physics and Astronomy, 
University of Manitoba, Winnipeg, MB, Canada R3T 2N2 



Abstract. We formulate the three-body problem in one dimension in terms of 
the (Faddeev-type) integral equation approach. As an application, we develop 
a spinless, one-dimensional (1-D) model that mimics three- nucleon dynam- 
ics in one dimension. Using simple two-body potentials that reproduce the 
deuteron binding, we obtain that the three-body system binds at about 7.5 
MeV. We then consider two types of residual pionic corrections in the dynam- 
ical equation; one related to the 27r-exchange three-body diagram, the other 
to the l7r-exchange three-body diagram. We find that the first contribution 
can produce an additional binding effect of about 0.9 MeV. The second term 
produces smaller binding effects, which are, however, dependent on the uncer- 
tainty in the off-shell extrapolation of the two-body t-matrix. This presents 
interesting analogies with what occurs in three dimensions. The paper also 
discusses the general three-particle quantum scattering problem, for motion 
restricted to the full line. 



1 Introduction 

The quantum mechanics of one-dimensional systems has always attracted interest 
in various research fields and has been extensively studied and investigated for 
a variety of purposes, ranging from the study of exactly solvable many-body 
systems |fl|, to the phenomenological description of tunneling phenomena @ 
in quantum wires, semiconductors, or across junctions H|. The mathematical 
properties for the underlying Schrodinger operators have been recently reviewed 
in ref. [Q] . Quantum scattering in one dimension has been also considered because 
of its value in pedagogical works ||, ||, 0, since it retains sufficient complexity to 
embody many of the physical concepts and behaviours that occur in the more 
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complex three-dimensional processes, but without the many technical aspects 
required for dealing with quantum scattering in higher dimensions. 

One-dimensional models are furthermore often constructed to gain deeper 
insight into the methods and approximations introduced to reduce the complexity 
of the higher-dimensional cases. In this spirit, we developed a one dimensional 
three-body model with the scope of testing three-nucleon theories in a simplified 
environment. We are in particular interested in the formulation of the quantum 
three-body theory on the full, (—00, +00) line, based on a system of coupled 
integral equations for the three-body transition operators (AGS equations |§, 
extended to an arbitrary number of bodies by Grassberger and Sandhas ||]). 
The approach is directly linked to the first mathematical theory of three-body 
scattering by Faddeev [pX| ], based on the integral formulation for the three-body 
Green's function and extended to N bodies by Yakubovski 11]. 

The formulation of the one-dimensional three-body problem in terms of the 
Faddeev- AGS scattering theory represents a somewhat "unusual" problem, since 
one-dimensional scattering corresponds also to the potential-barrier problem, 
which therefore has to be implemented with the clustering structure of the three- 
body problem. We illustrate in this paper an approach that harmonizes the 
typical two-channel structure of the "tunneling" problems, due to the presence of 
the transmission and reflection processes, with the cluster structure of the three- 
body problem, thus leading to an overall reduced S-matrix for the two-fragment 
processes which is 6 x 6. The formulation presented herein, with straightforward 
modifications, is suited also to describe the tunneling (or potential-barrier) effects 
due to the scattering of a two-body one-dimensional cluster impinging on an 
external barrier potential. A 3-D application along this line already has been 
discussed for a molecular system |L2| . 

The present work was motivated by the need to develop a simplified, 1-D test 
model for the non-relativistic treatment of the system consisting of three nucleons 
interacting via exchanges of a force-mediating "pion" (plus some shorter-range 
contributions). In modern three-nucleon dynamics, one of the crucial problems 
involves the treatment of those additional mesonic contributions that cannot 
be effectively described by the standard nucleon-nucleon potentials. Such terms 
are generally referred to as irreducible three-nucleon-force (3NF) contributions 
and are used to construct a phenomenological 3N potential. The hope is to cure 
in this way the discrepancies observed between experiments and rigorous Fad- 



deev three-body calculations with realistic 2N potentials |l^, 14]. The subsequent 
employment of the newer high-precision 2N potentials ]T5| , |16| , |17| ] essentially con- 
firmed the existence of such discrepancies. For example, the three-nucleon bound 
state, the triton, is typically under-bound by the high-precision NN potentials 
and the vector analyzing powers are substantially underestimated for low-energy 
nucleon-deuteron scattering (the so-called ^-puzzle) |18], p|. 

Considerable progress has been made since the introduction of 3NF based on 
isobars by Fujita and Miyazawa ||20|| , leading to a number of different 3NFs pi], 



22|, |23|, |2J] in use today, that can all be adjusted to describe adequately the 
triton binding energy. However, none of these terms have solved also the Ay- 
puzzle, which suggests that work still needs to be done to understand the role of 
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the 3NF terms in few-nucleon systems. A few years ago, an approach has been 
developed |^5| that reduces the field theoretic problem of the pion-three-nucleon 
dynamics to an approximate set of integral equations for the coupled 7rNNN- 
NNN system. The method extends the standard AGS formulation to the case of 
a three-body system with the additional pionic degree of freedom. The resulting 
equations merge together the Faddeev structure for the three-nucleon sector with 
the Yakubovski structure of the four-body sector, and has been shown to possess 
a kernel which is connected after iterations. 

Recently p6fl , a practical approximation scheme for the treatment of this 
coupled 7rNNN-NNN system has been proposed. The scheme is based on the 
assumption that the dominant part of the pion-exchange processes can be effec- 
tively described by contributions to the 2N potential, while the residual mesonic 
aspects lead to the addition of smaller, correction-type, 3NF effects. This assump- 
tion forms the basis for the standard description of nuclear systems in terms of 
aggregates of (pairwise) interacting nucleons, and ultimately relies on the chi- 
ral nature of the hadron interactions ^4|. The irreducible 3NF corrections due 
to the residual pion dynamics are treated perturbatively in the extended AGS 
equation, which is further reduced to an effective two-body integral equation in 
the inter-cluster momentum variable, by means of a finite-rank expansion of the 
two-nucleon potential. 

At the lowest order, the residual pion dynamics produced irreducible cor- 
rection contributions to the driving term of the AGS equation, and one type of 
these corrections can be identified with the so-called 2-7r-exchange 3NF diagram, 
and leads to the force models discussed, e.g., in refs. [20, 21, 22, 23, 24]. But 
at the same time, the approach produced also another type of irreducible pionic 



diagrams [gy, 27] - the l-7r-exchange 3NF term - that could not be identified 



with any previously known 3NF expression and which, as shown in completely 
realistic situations [28, 29], has the potential to solve the A y problem. For a brief 
review, see also ref. [30|. 

In the following, we develop a one-dimensional spinless model that mimics 
three-nucleon dynamics in one dimension. In particular, we investigate numeri- 
cally the effects of these irreducible "pionic" corrections in the one-dimensional 
three-body bound state by solving the homogeneous version of the AGS equa- 
tion. This 1-D model was originally constructed as a laboratory to investigate 
the dynamical effects of the meson on the three-body bound state [31]. However, 



one-dimensional systems in themselves are of considerable interest p2| , 33, 34 
and in many ways different from the three-dimensional case. This difference man- 
ifests itself, for example, in the existence of only two partial waves [@, H, HI and 



the different structure of Levinson's theorem |37], 38, p9| . 

Section |2| presents the AGS form of the three-body theory on the line, and its 
reduction to a Lovelace-type equation for a symmetrized three-particle system 
interacting via separable potentials. In Sect. ||| we define the explicit form of the 
pionic correction terms to the Lovelace-AGS equation. In Sect. |^ we consider 
the homogeneous version of the three-body integral equation, which leads to 
the calculation of the binding energies. In Sect. [5| we discuss the details of the 
three-body model and illustrate the results obtained. For the conclusions and a 



4 



Three-Body Dynamics in One Dimension: A Test Model 



summary of the results we refer to Sect. ||. 

Finally, we have also endowed this work with three appendices, devoted to 
specific aspects of the one-dimensional two-body system. Appendix A illustrates 
the integral formulation of the one-dimensional two-body scattering problem. Al- 
though the treatment can be found in the literature, e.g. in ref. [Jij, the discussion 
presented in the Appendix stresses the analogies with the three-body problem. 
Also the low-energy expansion of the two-body amplitudes is discussed therein. 
Appendix B illustrates the details for the expansion of the 2-body potentials in a 
finite-rank form, along the lines of the Unitary Pole Expansion (UPE) method. 
Appendix C discusses the method we have employed to solve numerically the 
one-dimensional two-body problem in one dimension, for both bound-state and 
scattering regimes. 

2 The Alt-Grassberger-Sandhas Three-Body Theory on the Full Line 

2.1 Distinguishable Particles 

In this section we formulate the Faddeev-AGS theory for three distinguishable 
particles moving in one dimension, (— oo, +oo). We consider three spinless parti- 
cles interacting only with 2-body potentials, and assume that these potentials are 
short-range, or at most exponentially decaying when the inter-particle distance 
goes to infinity. The case of irreducible 3NF-type diagrams will be discussed in 
Sect. ||. For simplicity, the two-body potential is assumed to support only one 
bound-state, as happens for the 2-nucleon system, but the formulation can be 
easily generalized to the case of a finite number of bound states. The particles are 
labeled "1", "2", and "3", and the indices a, b, c, ... range over the three labels. 
We will adhere to the odd- man-out notation which implies that when a = 1, for 
example, V a denotes the interaction between the pair (23). 
A system of three particles is described by the Hamiltonian 

H = H + J2v a . (1) 

a 

We then introduce the channel states, which are eigenstates of the channel Hamil- 
tonians 

H a = H + V a , (2) 
with a = 1,2, 3. The full Hamiltonian can also be written as 

H = H a + V a , (3) 

where 

V a = Y,V b (4) 

denotes the interaction "external" to channel a. In the one-dimensional system 
the channel states can be described by the following wave functions 



®±a = $±q a {Xa,ya) = e ±n * Va 4>a (x a ) , 



(5) 
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where 4> a (x a ) denotes the properly normalized bound-state wave function of the 
two-particle fragment and e iziqaVa denotes the relative motion of the particle a 
with respect to the cm. of the two-particle fragment (be). Here, x a ,y a are the 
Jacobi coordinates for partition a (be) and q a is the relative momentum between 
particle a and the cluster (be). The ± in the channel wave function is due to 
the presence of two physically different states, one with an incoming wave from 
the left and one with an incoming wave from the right. Note that this definition 
implicitly assumes q a ^ 0. Alternatively, we define also 

$a = $q a (Xa, Va) = e^<f> a (x a ) , (6) 

where q a may assume also negative values. 

The channel wave functions are eigenstates of the channel Hamiltonian H a 
with the eigenvalue 

h 2 q a 2 

where E a is the bound state energy of the two-particle fragment and M a is 
the reduced mass of particle a and the cluster (be). The masses of the three 
(identical) particles are given by m^c 2 = 938.27 MeV, the mass of the pion is 
m^c 2 = 134.98 MeV and the reduced mass of the particles b, c is [i a . Furthermore, 
because we consider three identical particles we also have \i a = [i b = \i c = \i and 
M a = M b = M c = M. 

Let us assume an asymptotic 1-D plane wave incoming from the left, de- 
scribed by the channel wave functions <?+ a - Then we can have three different 
transmitted two-cluster waves, and similarly for the reflected terms, which leads 
to the asymptotic form 

*Z (*, v) - h b6ab + R « b {E) ^ for Vb " "°° (8) 

where x, y is a generic set of Jacobi coordinates. The superscript L reflects the 
fact that the incident wave is coming from the left. Differently from the two-body 
case where there is only one single channel, we are now in a multichannel situ- 
ation where the rearrangement processes (i.e., a ^ b) have to be also described. 
This leads to the definition of the reduced S-matrix for two-fragment scattering 
which is formally equivalent to the one defined in Appendix A for two-particle 
scattering, but features in addition the cluster labels, 

Sah{E) -\R L ab (E) T*(E))- ^ 

This definition produces a 6 x 6 scattering matrix connecting all six possible "two- 
fragment asymptotic channels" #± a - The occurrence of such type of S-matrix has 
been observed already in the framework of an exactly soluble problem in ref. j5T] 
To solve the one-dimensional three-body scattering problem the S-matrix has 



to be determined, and for this we use the known expression [41] 



($a\ S ab |# 6 ) = 5 ab 5 (q a - q b ) - 2iri5 (E a - E b ) (<j> a \ V a \%) . (10) 
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It is then possible to introduce [^] the transition operators U ab , with the defi- 
nition 

<*a|^«6|*6> = ($a\V a \%) , (H) 

which satisfy the operatorial Alt-Grassberger-Sandhas equation 

U a b = SabGQ 1 + '^25 ac V c G c U c b , (12) 

c 

where S ab = 1 — 5 ab and G a (z) = [z — Ha)^ 1 . To recover the 6x6 structure for 
the reduced S-matrix we decompose the momentum on the line 

q = \q\ sign(g) , (13) 

where sign (q) denotes the sign of q. Then, we consider the two delta functions 
in the expression of the S-matrix, Eq. (|l~0|), 



Hqa ~ qb) = S (\q a \ - \q b \) <$aign( 9o ) sign(<fo) > ( 14 ) 

and 

5(E a -E b ) = ^8(\q a \ 2 -\q b \ 2 ) = JL-5 (\q a \ - \q b \) . (15) 



n01 J V\q b \ 

In this last equation we have assumed for simplicity that the bound-state energies 
of any two-particle fragments are the same. The generalization to the case with 
different binding energies, i.e. E a ^ E b for a ^ b, is straightforward. Therefore, 
the S-matrix can be written as 



(*«| S ab \$b) = 8 (\q a \ - \q b \) (*±«| S ab \$±b) , (16) 



where 



Sab (E) = ((f>± a \ Sab \4>±b) 

= x f 1 °\ g — /' t aft(+lgal. + lg&h- E ) tab(+\q a \,-\qb\;E)\ 
° ab [o i) V«6(-|?«|, + |<ft|;^0 *o6 " 1 ' 

To solve the AGS equation we expand the 2-body potentials into a finite-rank 
form. The method of expansion is discussed in Appendix B. For simplicity, in 
this section we illustrate the formulation for the case of a separable (rank-one) 
expression for the 2-body potentials 

sep _ Va\(t>a)((t>a\Va 

'\\Va\^a) ' [ ' 



where \<f> a ) corresponds to the properly normalized bound-state of the two-body 
subsystem a. Then the AGS equation transforms into an effective two-cluster 
Lovelace- type equation for the transition amplitudes [43], which is an integral 
equation in the inter-cluster momentum variable, 

Aab (<Za> q b \E) = Z ab (q a , q b ] E) 

/+oo 
Z ac (q a , q c ] E) A c (q c ; E) X cb (q c , q b ; E) dq c . (19) 
-oo 
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Figure 1. The exchange diagram, representing the AGS driving term 



The driving term of this equation corresponds to the exchange diagram repre- 
sented in Fig. [l], and is calculated according to the expression [44] 



Zab (Qa, %■, E) 



Xb 



2 J 



E 



2mjv 



<& + (Qa + qbf + (qb)' 



(20) 



Note that there is also a second driving term Z ac with a ^ c ^ b ^ a and 
implying the appropriate Jacobi-coordinate transformations. The A a (q a ;E) is 
given by 



A a (q a ;E) 



<t>a(p)Xa(p)dp 



Xa(p)Xa(p)dp 



oo E 



2M 



2fj, 



(21) 



In these last equations we have also introduced the form factors of the separable 
interaction, \\ a ) = Valfia), for convenience. 

Once evaluated for the on-shell momenta, the solution of this integral equa- 
tion provides the transition amplitudes 

{<t>±a\ U a b \<f>±b) = X ab {±\q a \,±\q b \;E) = t ab (±\q a \,±\q b \; E) , (22) 

which are needed in Eq. (|i~7|) for the construction of the S-matrix. 



2.2 Symmetrization 



Lovelace ]43| has discussed how the system of coupled equations can be sym- 
metrized, leading to a single-channel integral equation 

/+oo 
Z (q, q"; E) A (q"; E) X (q", q'; E) dq" . (23) 
-oo 

The q,q', and q" denote the inter-cluster momentum variables of the two- 
fragment states. The function A(q";E) is the same as calculated in Eq. (|2i"|), 
while the symmetrized driving term Z (q,q";E) is now twice that calculated in 



Eq. (20). The method works for both bosons and fermions, the difference be- 
tween the two being in the representation of the two-body sub-amplitude. For 
bosons, the two-body t-matrix has to be expanded via the UPE method using 
only symmetric states, i.e., with the even functions Xi(p) = Xi(~ p)- I n the case 
of three identical particles only one type of two-cluster partition is observable, 
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and therefore the 6x6 structure of the previously defined S-matrix and of the 
transition amplitude collapses back into a 2 x 2 structure, as was for the case of 
two-body scattering. 

The X (q, q'; E) represent the two-fragment transition amplitude which has 
to be evaluated, for a given energy E, at the only two possible on-shell values 
for q, q'. These momenta have to satisfy the condition 

9 ,2 2M . 2M - 

i = 4 2 = lf (E-E D ) = lf E, (24) 

where Ed is the bound state energy of the two-particle fragment. Because the 
system is defined on the real line, the values for the inter-cluster momenta have 
to be 



2M , . 2ME , . 

-w< E - E °)=i— ■ < 25 > 

This allows us to write the S-matrix for one-dimensional 2 + 1 scattering in the 
following way 

\0 1J h 2 \q\ \X(- \q\ , \q\ \E) X (- \q\ , - \q\ ; E) J 
which has the same structure as the two-body S-matrix defined in Appendix A. 

3 Treatment of the Irreducible Pionic Corrections 

The method, discussed up to now for a separable (i.e., rank one) expression for 
the two-body interaction, has to be extended to a rank-iV expansion, as described 
in Appendix B. The symmetrized Lovelace-type equation for a rank-iV potential 
is given by 

Xij (q,q';E) = Zy (q,q';E) 

N r+OO 

+ Y, Z ik {Q, E) A kl (q"; E) X tj (q", q'; E) dq" , (27) 

k,l J - 00 

where A^i is defined in Appendix B. The driving terms acquire a matrix structure 
in the rank space, Zij. 

Typically, the symmetrized Lovelace equation leads to a driving term 

Zij = M§, (28) 

which is twice that calculated in Eq. (^). The label "N" in the above expression 
is a reminder of the off-diagonal character of the driving term, due to the presence 
of 8 a b in Eq. (p!^). The factor 2 arises because in the symmetrized Lovelace three- 
body equation there are two off-diagonal elements that add coherently. 

The three-body Lovelace-type integral equation can be extended to include 



the degrees of freedom of one single meson [25]. To reduce the complexity of the 



resulting equation, one can consider a simplified formulation pq] where the effects 
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of the two-body interactions are treated exactly, and the residual irreducible 
(3NF-like) mesonic effects are treated approximately, as "first-order" corrections. 
The approximation scheme leads to mesonic corrections to the driving term, 
Eq. (p0|), of the standard three-body equation with two-body potentials discussed 
in the previous section. 

These pionic corrections corrections modify the driving term according to the 
expression 

Z ij = Z? + 2Z%, (29) 

where the superscript D denotes now the presence of diagonal-type contributions. 
Without these corrections one recovers the standard AGS-term by considering 
only the off-diagonal elements EI, i.e., Zf? = 0, and Zf, = Z, r,s 



In particular, the standard AGS term in the A-rank case is given by 
A ^,o , Xi (f + </) Xj \ -Q-i. 

Zg* (q, q', E) = I" V 2 > , (30) 

E-^^ + ^ + ^'f + ^'f 

where the form factors XiiXj are defined in Appendix B. We remind also that, 
since we have assumed a bosonic structure of the particles, only the even parity 
form factors are included in the separable expansion. 

The treatment of the irreducible pionic effects leads one to consider, in first 
approximation, the following correction terms 

zg = zt Gs + zir* 

Z§ = Zg"» (32) 

The diagrams generating these two additional terms are reported in Figs. ^ and||[ 
for z TPE3 and Z 3 , respectively (we recall that the standard AGS exchange 
diagram was reported in Fig. IT]). 




Figure 2. Correction terms to the off-diagonal AGS driving term 

The off-diagonal correction terms are given by 
Zjr™ (q,q>,E) = 

XdP) V™ 8 M,*«0 (33) 



v ' x p w (V) 

a 2M 2ju 2M 2/i 
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Figure 3. Correction terms to the diagonal AGS driving term. 



and similarly for the diagonal correction terms, 
Z% PE *{q, q ',E) = 

f - V OPE3 / / l\ X^) d d , (34) 

^ 2M 2/^ a 2M 2fj, 

The first contribution to the driving term Zij is generated by an intermediate 
7rN re-scattering process, as depicted by the diagram in Fig. ^. In the static 
approximation, the amplitude for such a process is given by the expression 

V {P,P,q,q) = — 7 ; "2^7 n : 72\ ' ( 35 ) 

Here, the Q, Q' are the pion momenta in the center of mass system with the 
proper set of Jacobi momenta 

Q=p-l-q> (36) 



Q' = ~q-P-^ (37) 

and i n N denotes the 7rN amplitude. 

One serious difficulty inherent to this test model was connected with the 
construction of this one-dimensional 7rN amplitude, since this represents a com- 
pletely abstract entity, which has no direct link to the phenomenological, three- 
dimensional amplitude of the real world []. We decided to circumvent this dif- 
ficulty by observing that the 7rN scattering lengths are roughly between one 
and two orders of magnitude smaller than the corresponding 2N scattering 
lengths Ea]. This suggests that probably the most simple parametrization for 
the 7rN amplitude in one dimension is to take the low-energy limit of the corre- 
sponding 2N one-dimensional amplitude, and to divide this amplitude by a factor 
ci between 10 and 100. The low-energy expansion for the two-body t-matrix is 

It should be noted that, in obtaining the 3NFs in the full three-dimensional form, this 7rN 
scattering amplitude is also the main source of uncertainty due to the need to find a satisfactory 
off-shell extension and suitable pion form factors |22|. 
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discussed in Appendix A, to which we refer for the details. The result leads to 
the following low-energy structure for the 1-D 7rN amplitude 

Un (Q, Q') ~ ~ [(a + &i) \Q\ \Q'\ + (& - &i) QQ'] , (38) 

where the parameter c\ has been arbitrarily set to c\ = —15, to be approximately 
consistent with the ratio between the observed scattering lengths of the 7rN and 
NN systems. 

Therefore, the final result for this contribution has the following structure 
{q,q,E) = -— (ao + ai) / - ■ — — ~ 2 _ — 2 j-^dp 



'j 



\Q'\ xjW) 



(q/)2 + (m*cy E _ wr _ (Ml 



2 

2M~ ~2fi 



f_ r _ , v f Q 

Cl7r «U J 2 (nhL£ ,2 ,. ,„„,-' ,•/„„- 



Xi {p) dp 



f Q' *M dv > (39 ) 

We stress that the arguments to arrive at the form Eq. (^) for the subtracted 
7rN amplitude have been ultimately quite arbitrary; however they reproduce the 
net result that, in the limit c\ — > oo, that is when the 7rN scattering length is much 
smaller than the 2N scattering length, the corresponding 3NF effect vanishes. 
The second contribution, v OPE3 , is generated by the diagram in Fig. [| The 
corresponding amplitude is obtained by summing over four terms, namely [^7]] 

V OPE *(p,p>, q , q >) = hG^G^ft + hG^G^ft 

+hG ( th 2 G ( i ) f{ + hG^Gfft , (40) 

where Gg 4 ^ is the 4-particle free propagator and /j denotes the creation operator 
of a pion on nucleon i, while /■ denotes the destruction of a pion on nucleon i. 
The quantity t\2 describes the two-body t-matrix for nucleons 1,2. Similarly to 
the case of the TPE3-term, also in this case one has to include a cancellation 
effect, since one must subtract from the two-body amplitude the corresponding 
"Born terms", i.e. the meson-exchange contributions, which can be identified 
with the two-body potential itself. Furthermore, the two-body t-matrix entering 
in this diagram has to be evaluated in a region which is substantially off-shell, 
because of the role of the exchanged meson in this diagram. Therefore, following 
the studies of refs. 27, 25, 2£] we have considered the effective parameter that 
governs the cancellation for the diagonal correction terms: 

ti2 {p,p; z) = c 2 t 12 [p,p'; z) - v u (p,p') , (41) 

where 

z = E-^L-?uuAQ) ■ (42) 
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The effective parameter C2 represents an overall correction factor for the far- 



off-the-energy-shell extrapolation of the t-matrix [28, 29 1. Also, V12 denotes the 
meson exchange contributions as defined in the nucleon-nucleon potential acting 
between nucleons 1,2 and 



<rrj2 s>2 

uAQ) = \jQ 2 + ^- (43) 

is the relativistic pion energy with the pion momentum 

Q = q-q'. (44) 

Using a static approximation for the four-body Green's functions and in- 
serting the explicit expressions for the creation and destruction operators gives 
the following form for the static contribution to the diagonal correction terms 

v OPE3 (p, P ', q , q >) fl 2 



it hcuj^ (Q) 

1x2 (p,p; E - - hcu* (Q)J + ti 2 Up'; E - - hao* (Q) 



(45) 



In this paper we are using the unitary-pole expansion and consequently the 
subtracted NN t-matrix is defined as 

t 12 [p,p'\ E - M. - hcuj n (Q)j = 



Y, Xi (p) ( c* A? ( E ~ w - ^ (Q) ) + %% I Xi [p') , (46) 



2 



2M 



which leads to the final expression for the diagonal contribution of the correction 
term in a rank-iV separable expansion 



7 OPE3( i jp\ f 2 2 I ^ f Xi (p) Xk (p) j 

P2^« ( # " ^ " (Q) J + c 2 A kl (e - - haon (Q) j + 2^^ 



x 



xKgWgO (i. 7) 



2M 2fi 



This fully defines the driving term of the rank-N symmetrized Lovelace-type 
equation. 
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4 Bound-state equation 

In the previous sections we have provided the AGS formulation for three-body 
scattering in one dimension. However, we will now restrict the discussion to 
the corresponding bound-state problem. We will apply this formalism to the 
calculation of the binding energies, as a first test for analyzing the effects of 
the irreducible pionic corrections in the 1-D three-body system. As is known, 
the Faddeev or AGS method formulates the scattering problem in terms of an 
inhomogeneous integral equation which can be reduced, after iteration, to one 
of Fredholm type. Such equations are known Era to have no solution in case the 
related homogeneous equation admits non-trivial solutions, which are associated 
with bound states. 

A further step, introduced for the convenience of numerical solution of the 



homogeneous AGS-Lovelace Eq. (27), is to express it in the form of a sturmian 



eigenvalue problem [[47]]. Therefore, in order to find the three-body bound state 
energy, we solve the generalized eigenvalue equation [48] 

J Utj (q, q', E) i> jk (q') dq' = Vk (E) ^ J Z {j (q, q' , E) (q') dq' . (48) 

3 3 

In this equation 



Uij 



(q, q', E)=J2 f Z ik [q, <l", E) A M (q"; E) Z Xj (q", q' , E) dq" , (49) 
kl 



and the remaining quantities are the same as those defined in Sect. |3| The eigen- 
values r\ k (E) are energy-dependent, because of the dependence on the energy 
of Z and A. The energy is varied, and the bound-state energies of the three- 
body system are found by the condition that the eigenvalue be equal to one at 
E = Eb- Since the binding energy is negative, all integrals (in Sect. ||) involv- 
ing the Green's function are nonsingular, while for energies above the scattering 
threshold, whenever such integrals become singular, they have to be regularized 
with suitable subtraction techniques. 



5 Results of the Test Model 



We have constructed three different two-body potentials that mimic the nuclear 
force on the line. The general structure of these potentials is given by 

V (x) = V A <T m + V R e~ n ^ x \ , (50) 

with the free parameters Va, Vr. In Table [l] we report the three different sets of 
parameters. The range parameter (3 of the attractive part is the inverse of the 
Compton wavelength of the pion 

0=^1 (51) 

which we consider fixed. The range parameter of the repulsive part is variable and 
has been parametrized in terms of multiples of (3, denoted as n. That is, n denotes 



14 



Three-Body Dynamics in One Dimension: A Test Model 



the ratio between the Compton wavelength of the "scalar pion" , the boson that 
mediates the attraction component, and that of the repulsive short-range com- 
ponent. For the purpose of the upcoming discussion it is unimportant whether 
the short-range part is mediated by heavy-boson exchange, or by a generic cut- 
off mechanism which suppresses the contributions coming from momenta higher 
than the cut-off parameter. Three values for n have been selected, 5, 10, and 20, 
corresponding to an increasing "stiffness" of the short-range repulsion. 

The remaining parameters were fitted to reproduce a two-body bound state 
at (approximately) the deuteron energy and a two-body low-energy scattering 
parameter at about 24 fm, which would then correspond to the 1-D equivalent 
of the NN scattering length. The binding energies reported in Table |l| were 



Table 1. Strength and range parameters of the one dimensional poten- 
tials approximating nuclear physics and the imaginary parts of the low 
energy expansion parameters 





n 


Va [MeV] 


Va [MeV] 


do [fm] 


di [fm] 


E B [MeV] 


v 5 


5 


-24.51 


65.717 


24.46 


-14.22 


-2.225 


Vio 


10 


-23.71 


129.017 


24.43 


-14.18 


-2.225 


V 20 


20 


-23.49 


262.321 


24.40 


-14.17 


-2.225 



calculated using the coordinate-space method described in Appendix C. The low- 
energy scattering parameters &o,&i, introduced in Appendix A, were calculated 
numerically with the same method. 

These potentials have a simple representation also in momentum-space 

V (k, k>) = ?-( % + ,) (52) 

V ; vr \(k- k'f + (3 2 (k - k'f + {n(3f J K ' 

and it has been checked that the numerical calculations done in momentum space 
are consistent with the results obtained with the coordinate-space method. 

The bound-state energies of the three-body system are found by the condition 
that the generalized eigenvalue in Eq. (p8|) has to be equal to one. We investigated 
the effect of the correction terms for the three potentials V5, Vio, and V2Q. Note 
that the coupling constant of the force-mediating boson is implicitly defined 
by the two-body potential, since the boson-exchange potential is given by the 
exponential 

f je-^. (53) 

This is the analogue of the Yukawa potential in one dimension. Hence the cor- 
rection terms discussed in section |3] have to be evaluated with the condition 
f = (3V A - 

In Fig. I we show the convergence behaviour of the three-body binding energy 
in respect to the rank of the separable expansion for the three potentials. The 
left graph shows the convergence of the calculation without any correction terms 
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2 3 4 5 6 
Rank 



Figure 4. Convergence behavior of the binding energy for the three potentials in respect to 
the rank of the UPE. The left graph shows the calculation with no correction terms and the 
right one with both correction terms (diagonal and off-diagonal) included. The o denote the 
results for potential Vs, the x for potential Vio and the + for potential V20. 



and the right one with both corrections (diagonal and off-diagonal). The two 
effective parameters c\ and c<i have been set to —15 and 1, respectively. It is 
clearly seen that both calculations show a convergence in respect to the rank 
of the UPE. In Table ^ we show the numerical results for rank 7 calculations 

Table 2. 'Triton' energies in MeV for the three potentials, including 
different 3NF corrections. The results are given for a rank 7 UPE and 



with 


an effective parameter C2 


= 1 






2NF 


2NF+OPE3 


2NF+TPE3 


2NF+OPE3+TPE3 


V5 


-7.722 


-7.708 


-8.733 


-8.718 


Vio 


-7.684 


-7.662 


-8.670 


-8.650 


V20 


-7.668 


-7.642 


-8.631 


-8.608 



with different combinations of correction terms. It is worth noting that, once the 
local two-body potential has been constructed to mimic the properties of the 
three-dimensional 2N system (2N binding energy and low energy parameters), 
the three-body binding energy turns out remarkably close to the values expected 
for the 3N system, in spite of the fact that we have considered a rather different 
system, namely three spinless bosons defined in one dimension. 

It has also to be observed that the off-diagonal pionic contributions, orig- 
inated by the TPE3 diagram of Fig. |l| can indeed be modelled to correct a 
hypothetical under-binding of the three-body system, since these terms have 
the potential to provide an additional contribution to the binding energy which 
can be sized around 0.9 MeV. This also presents interesting analogies with the 
3-D case of realistic 3N system. We stress, however, that details can change, 
also significantly, depending on the particular model considered for the off-shell 
extrapolation of the 7rN amplitude. 

The diagonal correction terms related to the OPE3 diagram (Fig. [|) produce 
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changes in the three-body binding energies which are significantly smaller than 
the off-diagonal ones. This also presents interesting analogies with the realistic 
3N case [ 28 [ , where it has been already found that the triton binding energy is 
not strongly affected by the diagonal correction terms. Within this 1-D model 
one could conclude that they could be neglected in first instance, if the effective 
parameter c<i is set to one. 

Of course, this should not lead to the assumption that the OPE3 terms can 
also be ignored in calculations of other observables, as it happens for the case 
of the vector analyzing powers, for realistic nucleon-deuteron scattering at low 
energy. It was shown in refs. [28, that the effective parameter C2 can be fitted 
to the peak of the Ay-data. The value of the parameter in these fits depends 
sensitively on the particular NN-potential used in the calculations, but in every 
case the resulting curves reproduce the existing data base quite well. Finally, it 
is worth mentioning that the two correction terms produce additive effects to 
the three-body binding. 



0.8 1.0 1.2 

effective parameter c 2 



Figure 5. Dependence of the binding energy on the parameter ci in a rank 1 calculation for all 
three potentials including the diagonal correction terms. The o denote the results for potential 
Vs, the x for potential Vio and the + for potential V^o- The thick line denotes the triton energy 
without any correction terms for all three potentials. 



In Fig. U we show how varying the parameter C2 affects the binding energy 
in the one-dimensional three-body model. The range of variation for C2 shown in 
the figure is consistent with the values for C2 found in refs. [28, 29 1. As mentioned 
earlier, one possible reason for using C2 as an adjustable parameter is the fact that 
the 2N t-matrix entering in the OPE3 diagram has to be calculated off the energy 
shell, in a region which is difficult to constrain by experimental data. In addition, 
the figure reveals how the variations with C2 are dependent also on the selection of 
the two-body potential. At a parameter value of 0.7 the diagonal correction terms 
introduce an effect that is approximately 10% of the effect introduced by the off- 
diagonal correction terms. At a value of 1.3 the effect is in the opposite direction 
and the magnitude of the effect is strongly dependent on the particular potential, 
as seen by the different slopes. We have chosen to show the rank 1 calculation for 
different potentials, because the binding energies without correction terms are 
practically equal for all three potentials in this case. However, the same behavior 
can also be observed in the higher-rank calculations. We conclude that, a priori, 
one should not neglect the diagonal terms in this model calculation. 
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6 Summary and Conclusion 

The scope of this study was originally the development of a strongly simplified toy 
model in which to test a three- nucleon theory with residual pionic corrections, 
namely with contributions that could not be described via conventional two- 
nucleon potentials. The general theoretical framework adopted as starting point 
has been the coupled equations of ref. p3] for the 7rNNN-NNN system, further 



reduced in ref. [26] to an approximate but more practical three-body dynamical 
equation with irreducible pionic corrections. 

For this scope we formulated the quantum mechanical three-body problem 
in one-dimension, on the (— oo,+oo) line, and considered the addition of two 
different 3NF-type diagrams generated by the residual pion dynamics. In the 
specialized literature, attention has been focused - also recently - to these types 
of diagrams in realistic 3N studies for both bound and scattering regimes. 

In this work, we concentrated our numerical study on the binding energies of 
the three-body system. We described the gross features of the 3-body bindings, 
observing also interesting analogies with the more realistic 3N-bound states, de- 
spite the rather severe simplifications introduced by the one-dimensional model. 

In particular, using simple 1-D potentials that could reproduce the deuteron 
energy, we obtained a 3-body binding of about 7.5 MeV. Addition of residual 
pionic diagrams of the "TPE3" type (Fig. |2|) carries an additional binding of 
about 0.9 MeV. However, significant model dependence on the structure of the 
"7rN" amplitude (related to these diagrams) can and have to be expected. The 
addition of the second type of diagrams, denoted "OPE3" (see Fig. |||), produces 
an additional effect on the binding which is either negligible or moderate, de- 
pending on the specifics of the selected two-body potential, and on the value of 
an effective parameter that governs such diagrams. It must be observed that one 
cannot conclude that the OPE3 terms are less important, since they could quite 
significantly affect other three-body observables. This indeed turned out to be 
the case for the nucleon-deuteron A y observable in realistic studies [28, 29 1. 



However, we realized also that the 1-D approach we developed for the scopes 
discussed above went somehow beyond the initial intentions of the authors. In 
fact it turned out, to the best of our knowledge, that the three-body quan- 
tum scattering theory using the Faddeev-AGS-Lovelace formulation in terms of 
coupled integral equations has never been discussed before for one-dimensional 
systems, while the general formulation for such systems could attract interest 
in various fields dealing with one-dimensional systems. For this reason, we have 
devoted Sect. |2| to the general discussion of the quantum-mechanical 3-body scat- 
tering problem in one dimension. Both cases involving distinguishable or identical 
particles have been formulated in terms of the integral-equation approach. 

As possible applications and future extensions of the 1-D model considered 
herein, one could consider also other cases of three-body models where there 
is the need of an explicit coupling to the meson degree of freedom. This is the 
case, for example, of the constituent quark model based on Goldstone-boson 
exchanges [49], which provided a rather realistic description of the excitation 



spectra of all light and strange baryons [50]. One might think to improve the 



description of the excitation spectra by including explicit mesonic degrees of 
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freedom into the resonance states, and in this framework to study the tt and 
77 decay modes of N and A resonances. Another possible application is in r/- 



3N scattering [51] and associated reactions [52|. These processes are dynamically 
coupled with the absorption channel, and with the 7r-channel itself, and therefore 
the dynamics is highly complicated. In both these examples, it would be rather 
useful to develop similar 1-D models as theoretical laboratory where to test 
methods, approximations, and dynamical assumptions in a simplified framework. 
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Appendix A: Partial Waves and Low Energy Behaviour in One- 
Dimensional Scattering 

Here we review some general results obtained in 1-D scattering theory. We consider a particle of 
mass m with no internal degrees of freedom subject to a potential V (x) in one dimension. The 
potential is assumed to be real and approaching zero sufficiently fast for x — > ±00. Furthermore, 
we restrict our investigation to symmetric potentials V (x) — V (—x). The Schrodinger equation 
for this problem is 



dx 2 

where k = (2m_E) 5 /h is the wave number, E is the energy of the particle and U (x) — 
(2m/fi 2 ) V (x). In this section we examine the transmission-reflection problem for positive en- 
ergies. This problem has two independent solutions with incident waves from the left and right 
respectively. The asymptotic form for the wave function incident from the left is given by 

ViM^l^+toWe-*-, for X ^ -00 

V ' [T L (E)e lkx , for x — > +00 . K ' 

The transmission and reflection coefficients Tl,Rl lead to the transmission and reflection prob- 
abilities, which have to satisfy 

\T L {E)\ 2 + \R L (E)\ 2 = 1 (A.3) 

due to conservation of probability. The asymptotic form for the wave function incident from 

the right is defined similarly replacing x with —x. 

The S-matrix is given according to the definition j3f], ^] 

c,™_(Tl{E) Rr{E)\ 

which reduces to the identity matrix when the potential goes to zero. In general, the matrix 
elements are complex and hence the matrix is defined by eight real parameters. Unitarity of the 
S-matrix leads to four constraints on the matrix elements and for real and symmetric potentials 
two more constraints are observed. This leads to the identity of the transmission and reflection 
coefficients for the incident waves from the left and right 

T (E) = T L (E) = T R (E) (A.5) 

R (E) = Rl (E) = Rr (E) , (A.6) 
which represent the only two independent parameters of the S-matrix 



c (F) -{T(E) R(E) 
1 ' \R(E) T(E) 



(A.7) 
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In the one-dimensional scattering problem there are two partial waves, one with even parity 
and one with odd parity. The S-matrix in respect to these partial waves for real and symmetric 
potentials is found by a diagonalization procedure 



S(E) = 



T(E) + R (E) 






T(E) -R{E) 



(A.8) 



Furthermore, due to the unitarity of the S-matrix, it can also be given in the partial wave 
representation 



S{E) = 







(A.9) 



where r5n,i {E) are real and define the phase shifts of the scattering process [^, ^4j. In Fig. A..1 
we show the odd and even phase shifts calculated with the procedure described in Appendix C. 
They show the behaviour expected from the one-dimensional Levinson theorem [B7L p8| , namely 



0.001 0.01 0.1 1 10 100 1000 
Energy [MeV] 



Figure A.l. Phase shifts for the odd and even wave function with potential V20, + denotes 
the odd phase shift and x the even phase shift. 



the odd phase goes to zero at threshold and the even to ^ , because of the even bound state. For 
completeness, we note that the transmission and reflection coefficients can also be expressed in 
respect to the phase shifts in the following way 

T(E) = l( e ^oW +e HiW) = cos (^o (E) - S 1 {E)) e^ {E)+ ^ E)) (A.10) 

R(E) = i( e 2i5o(B) -e 2Wl(B) ) =sin(5 (i5)-<5i(i5))e i ^ (E)+ * l(B)+ f). (A.ll) 

Finally, according to Galindo and Pascual |J, for potentials that go to zero sufficiently fast, 
when x goes to plus or minus infinity, the following generic behavior is observed 



*<*> » W 2 ?. {^-» < A - I2 > 



R(E) » -l + Qi y — , ^—^0, (A.13) 

where an, cei are generally complex constants not equal to zero. Therefore, the S-matrix in the 
low-energy limit tends to a total reflection. It should be noted that, for the specific type of 
potentials we use, it turns out the constants <an,ai are purely imaginary numbers. 

The transmission and reflection coefficients c an be calculated by resorting to the Green's 
function formalism |l[ The Schrodinger Eq. ( [A~l| ) can be written in integral form, which 
leads to the Lippmarm-Schwinger equation 

/+00 
G (x,x')U(x')V k (x')dx' , (A.14) 
-OO 
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where $ k (x) is a solution of the free Schrodinger equation and the Green's function Go (a;, x') 
is defined incorporating the boundary conditions of the equivalent Schrodinger problem. For 
1-D systems the Green's function with outgoing boundary conditions has the following form 



Go {x,x') = ~^^ X ~ X '!. 
The general solution of the free Schrodinger Eq. ( |A.l| ) is given by 

(x) = Ae ihx + Be~ ikx , 



(A.15) 



(A.16) 



where the coefficients A and B are arbitrary. Inserting these expressions in the Lippmann- 
Schwinger equation leads directly to the following equation 



&h (x) 



Ae kx + Be 1 x 



im 



r e ik( . x - x ")V(x')V k (x')dx'+ f + °° e ik ( x '~ x )v(x')V k (x') dx' 



(A.17) 



and A, B can now be determined according to the various boundary conditions. One choice is 
A = 1, B = and corresponds to the following equation 



<{x) = 



h 2 k 



im 



ikx 



-ikx 



V (x) (x') dx' 

V (x) &k {x') dx' 



(A.18) 
(A.19) 



which is equivalent to the previously defined solution in Eq. (A.2). The transmission and re- 
flection coefficients for scattering with an incoming wave from the left are then given by the 
expressions 



where 



rkL 

J± 



Tl(E) = l + ^ft L 

Rl(e) = ±f!L\ 



2m 



^1 e^ xk V(x')$£(x')dx' 



(A.20) 
(A.21) 

(A.22) 



is the scattering amplitude for an incoming wave from the left. In a similar way it is possible to 
describe the scattering system with an incoming wave from the right by choosing the coefficients 
A = 0, B = 1. However, we already argued that for the type of potentials we use in this paper 
the left- and right- coefficients are equal. Therefore, we drop the indices L, R in the remainder 
of this section and denote the scattering amplitude in both directions by f± . Observing that in 
the one- dimensional system, for a given energy E, there exist exactly two possible momenta 

k±=± V2mE 



we define the on-shell t-matrices by the expression 

t(k ± ,k;E) = (k±\t(E)\$ k ) = {$ ±k \V\* k ) 

e^ kx V{x')* k {x')dx' = -^f k ± . 

For the off-the-energy-shell situation the t-matrix is introduced according to 

t(E)\<P q ) = V\<P q ) 
and it can be calculated by the Lippmann-Schwinger equation 

t (E) = V + VG t (E) 



(A.24) 

(A.25) 
(A.26) 
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which reads in momentum-space 

t(q,q';E) =V(q,q') + 
where 

Go (q,E) 



V(q,q") G (q",E)t(q",q';E) l 



E 



j-2 2 \ -1 

n q x 
2M 



(A.27) 



(A.28) 



Finally, combining Eqs. ( |a~^ , |A.2C| , [a!2i| , |A.24| ) it is possible to express the S-matrix in respect 
to the on-shell t-matrix 



S(E) 



I 


; 


im 1 






~Wk\ 



im ( t(k,k;E) t(k,-k;E) 



(A.29) 



A diagonalization procedure lets us express this result in terms of two "partial wave" t-matrices 



S(E) 







im / 


G 




~ hPk \ 



where 



to (E) 
ti (E) 



■(E) 
ti (E) 



t(k,k;E) + t(k,-k;E) 
t(k,k;E) -t(k,-k;E) 



(A.30) 



(A.31) 



The low energy behavior of the t-matrix can be defined by a direct comparison to the low 
energy behavior of the on-shell S-matrix, which yields 



to [E) 



ti (E) 



im 
im 



2mE 

h 2 

2mE 



(a + on) 



(ao — Qi) 



2mE 



2mE 



(A.32) 



(A.33) 



where 0, 1 denote the even and odd t-matrix respectively. In Fig. A. 2 we show the low energy 
behaviour of the transmission and reflection coefficients explicitly. The figures were calculated 
using the procedure outlined in Appendix C. Clearly, the imaginary parts of the reflection and 



' - XXX* vyxX>X' XXXVX x. . . .v ;7 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 
k[fm" 1 ] x10 -3 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 
k[fm" 1 ] x10" 3 



Figure A. 2. Real and imaginary parts of the coefficients 72+1 and T . The + denote the 
imaginary parts of the coefficients, the x denote the real parts of the coefficients. 



transmission coefficients show the linear behaviour, while the real parts are equal to zero. In 
the TPE3 a subtracted 7rN t-matrix enters and we approximate it by the definition 



tF (q,q) = 



(ao + «i) |g| \q'\ + (ao — «i) qq' 



(A.34) 



Here, the param eters &o,&i denote the imaginary parts of ao,ai, wh ich are given by the 
Eqs. (A. 12 A. 13). Furthermore, to get this result from Eqs. (A.32 A.33), it must be observed 
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that both odd and even terms must be added coherently, since the meson and the nucleon are 
distinguishable. Fin ally, w e have neglected the additional contribution due to the presence of 



the term —1 in Eq. (A. 13). This contribution originates because in the low-energy limit the 1-D 
scattering event becomes a purely reflective process (T — > 0, R — > —1). It is eliminated from 
Eq. (^) because in the 7rN t-matrix there are contributions that need to be subtracted, i.e. the 
forward propagating nucleon diagram and its crossed counterpart. In the three-dimensional case, 
the subtraction of such diagrams is discussed for example in ref . . Finally, a phenomenological 
constant Ci has to be introduced to account for the different magnitudes of the low energy 
7rN and NN t-matrices. For insertion into the off-diagonal correction terms the expression is 
rewritten in the appropriate units multiplying it by the constant (2m) /fi 2 . 



Appendix B: The unitary pole expansion 

In this appendix we describe the Unitary Pole Expansion (UPE), which is the input for our 
three-body bound state problem. The UPE is closely related to the Weinberg series and is 
described by Harms pffl. We derive the UPE in a slightly different, but equivalent, way as dis - 



cussed by Pisent et al7[B8|. The homogeneous form of the Lippmann-Schwinger Equation (A. 26) 
can be given by the generalized eigenvalue equation 

VG (E) V \<f> n (E)) = Vn (E) V \<f> n (E)) . (B.l) 

The form factors of the separable expansion are defined by 

\ X u(E)) =V\<t> n {E)) . (B.2) 

In our case we have only one bound state and we choose to use the expansion at the bound 
state energy E = — B. In this way we recover the Unitary Pole Expansion (UPE). Furthermore, 
it is easily seen that the first eigenvalue 771 (E) is equal to one in this case. Therefore, the bound 
state \4>b) and the first form factor in the expansion satisfy the relation 

\xi)=V\4>b). (B.3) 

The bound state calculated with this definition was compared to the one calculated with the 
Fox-Goodwin method and it was seen t hat they indeed do agree with each other. It is also 



known [ft8| that the eigenvectors of Eq. (B.l) have the following orthogonality relation 



i\ V\<t> m ) = -Vn$nm , (B.4) 



where the \</> m ) are the solutions of Eq. (B.l ) at the binding energy E = —B. This orthogonality 
relation also specifies the normalization of the eigenfunctions, which is slightly different from 
that introduced in Section |^. The two-body potential V is given in respect to the following 
separable expansion 

V = -Y] lx ' )<Xt| (B.5) 

and the separable T matrix is given by 

r = X)|xi)4*i(Xil, (B.6) 

where the matrix elements of A are defined by the expression 

[A-\ = -TH6i i -(xi\Go\Xi) • (B.7) 

In the case that only one term in this expansion is retained the expansion reduces to the Unitary 
Pole Approximation (UPA) discussed by Lovelace ^3] and Fuda [^| 

TupA = \xi)AutXi\ (B.8) 

with 

^11 = -[l + <Xi|G |xi)] _1 • (B.9) 

It is easily seen that both the UPE and UPA have the right pole structure, namely the separable 
T matrix has a pole at the bound state energy. 
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Appendix C: Numerical Studies of the Two-Body Problem in Coor- 
dinate Space 

We note that our one-dimensional system is defined on the whole line and not the half line like 
the radial Schrodinger equation of a three-dimensional problem. In this Appendix we show how 
the system can be solved numerically for any finite-range symmetric potential. 
The problem is defined by the time-independent Schrodinger equation 



which can be written as 



where 



HtP (jc) = EiJj (x) , (C.l) 
u" (x) + w (as) u (x) = , (C.2) 



w(x) = ^(E-V(x)) . (C.3) 



The Fox-Goodwin method j58| replaces the differential equation by a finite expression 

2u (x) — u(x — h) — (jz) (Ww (x) u(x) +w(x — h)u(x — ft)) 
u(x + h) = V / / ,, N , (C.4) 



2u (x) — u (x — ft) — 




(lOw (x) u (x) +w(x — h)u(x— ft)) 


1 + 


I 18; 


w (x + ft) 



which is exact to the order of ft 4 . If one knows the values of the function u at the positions 
and ft, then it is possible to calculate the value of the function at the point 2ft. Repeating this 
procedure with the function at ft and 2ft will give its value at 3ft and so on. With this procedure 
the function is defined on a grid with discretization interval ft. The first two values chosen in 
this procedure, namely u (0) and u (ft) depend on the boundary conditions of the system. This 
procedure works well for the radial equation, because one has a definite starting point at zero. 

The generalization to the full line is simplified if one takes advantage of the symmetric 
nature of the potential. The symmetry implies that the eigenstates can be classified into sym- 
metric (even) or anti-symmetric (odd) states. Therefore, it is possible to restrict the calculation 
to the positive part of the real line by imposing the appropriate boundary conditions at the 
origin for both odd and even states. 

Odd states have to go through zero at the origin. Consequently, we define the boundary 
conditions at the origin in the following way 

Mi (0) =0 wi (ft) = 1 . (C.5) 

Even states have to go through an extremum at zero, but without crossing the origin. Conse- 
quently, we define the symmetric boundary conditions as follows (h — > 0) 

u (0) = 1 u (ft) = 1 . (C.6) 
C.l The Bound State Problem on the Real Line 

Bound-state solutions, in both odd and even cases, have to be normalizable. In the positive 
direction, negative-energy solutions follow the asymptotic behaviour 

u (x) = ae~ KX + be KX , (C.7) 



with k — sj2m \E\/h being the wave number for E < 0. 

The normalization condition determines the bound-state by selection of the energies cor- 
responding to a zero of the b coefficient, namely a vanishing wavefunction for x — > +oo. Nu- 
merically, the bound states and energies are fo und by s tarting the Fox-Goodwin algorithm at 



x — with either the odd and even Conditions (C.6 C.5| ), and then evaluating the wavefunction 



for very large x valu es. G enerally, the function assumes very large values, due to the increas- 



ing exponent in Eq. (C.7), except for energy values corresponding to bound states, where the 



wavefunction drops to zero at large distances and becomes normalizable. 
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C.2 The Scattering Problem on the Real Line 

The 1-D scattering problem is solved separ ately for ant isymmetric and symmetric states by 
starting at the origin with either Condition ( 3.t ) or ( |C.5|) ) and propagating the solution in the 
asymptotic region x — > +00 with the Fox-Goodwin method. In the odd case we know that the 
solution in the asymptotic region has to be of the form 

wi (a;) = ai sin (kx + Si) (C8) 

with k = \/2mE/h being the wave number and the bar denoting the asymptotic solution. 

Equating the last two points of the numerical solution with this general expression gives 
the following two equations 

ui (B) = ai [cos (Si) sin (kB) + sin (Si) cos (kB)] 
Ui{B-h) = ai [cos (Si) sin (k (B — h)) + sin (Si) cos (k (B — h))] , (C.9) 

where B and B — h are the largest and second largest :r-value at which the wavefunction is 
calculated numerically. Solution of this system leads to the following expression defining the 
phase shift 

tnn ( g ) _ sin [ fc ( B ~ ( B ) - sin i kB } "1 (g - h ) / rim 

1 1 ' — cos [k (B — h)] ui (B) + cos [kB] ui (B - h) ' 1 ' 



Substituting the phase shift back into (C.£) gives also the odd amplitude ai. 

The even phase shift can be similarly determined by the asymptotic behavior 

wo (x) = ao cos (kx + So) , (C.ll) 

which yields 

(x \ — cos \kB] up (B - h) - cos [k (B - h)] up (B) . . 

anl 0j ~ sm[kB]u (B-h)-sm[k(B-h)]u Q (B) ' 1 ' 

As for the antisymmetric case, also the symmetric amplitude ao can be determined in the same 
way. 

Generally, the asymptotic behavior of the one dimensional scattering system for incoming 
waves from the left (x — > — 00), is described by 

u{x^-oo) = e ikx + R ■ e~ ikx (C.13) 
u{x^+oo) = T-e lkx , (C.14) 

where R and T denote the reflection and transmission coefficients. We define the scattering 
wavefunction u(x) with the proper asymptotic behaviour by a (suitably chosen) linear combi- 
nation of odd and even solutions, namely 

u (x) = up (x) + i ui (x) . (C.15) 

ao ai 



Comparing the asymptoti c behaviour for both antisymmetric and symmetric Solutions C.8 C.ll 



with the Form C.13 



C.14 one derives the two scattering coefficients 



T(E) = i(e 2iWE) +e 2i<lW ) (C.16) 

R(E) = - ( e ™ W_ e 2iS 1 (E)\ (C17) 



2 

according to the expressions already mentioned in Appendix A 
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